- 9-3 


NASA 
technical 
iPaper " 
3350 


fJ "> y /C 


K. 3 


August 1993 


Semiempirical Fragmentation 
Models on Galactic Cosmic 
Ray Transport Calculations 
With Hydrogen Target 


Judy L. Shinn, 

John W. Wilson, 
and Francis F. Badavi 


(NAS A- TP-3350) SEMIEMPIRICAL 
FRAGMENTATION MODELS ON GALACTIC 
COSMIC RAY TRANSPORT CALCULATIONS 
WITH HYDROGEN TARGET (NASA) 18 p 


N94-13263 


Unci as 


Hl/93 0182A94 







nil M IIH1II Mil nil II lil liaa 




NASA 

Technical 

Paper 

3350 


1993 


Semiempirical Fragmentation 
Models on Galactic Cosmic 
Ray Transport Calculations 
With Hydrogen Target 


Judy L. Shinn 
and John W. Wilson 
Langley Research Center 
Hampton, Virginia 

Francis F. Badavi 
Christopher Newport University 
Newport News, Virginia 


NASA 

National Aeronautics and 
Space Administration 

Office of Management 

Scientific and Technical 
Information Program 



Abstract 


Nuclear fragmentation cross sections of Silberberg and Tsao that 
are more accurate for a hydrogen target have been implemented in 
the data base to replace those of Rudstam for a galactic cosmic ray 
transport code (HZETRN). Sample calculations have been made for the 
transported galactic cosmic ray flux through a liquid hydrogen shield 
at solar minimum condition to determine the effect of such a change . 
The transported flux based on the Silberberg- Tsao semiempirical for- 
malism contains fewer high-LET (linear energy transfer) components 
but more low-LET components than the results based on Rudstam s for- 
malism,; and this disparity deepens as the shield thickness increases. A 
comparison of the results obtained from using both energy- dependent 
and energy-independent cross sections of Silberberg and Tsao indicates 
that the energy- independent assumption results in an underestimation 
of high-LET flux above 100 keV/fim by approximately fO percent for 
a 15-g/cm 2 thickness of liquid hydrogen . Similar results were ob- 
tained in a previous study (NASA TP-3243) when both energy- dependent 
and energy-independent cross sections of Rudstam were considered. 
Nonetheless , the present study found that an energy-independent cal- 
culation would be best accomplished by using Rudstam f s cross sections 
as done in the past for various engineering applications. 


Introduction 

Estimates of galactic cosmic ray exposure and 
the required shielding for astronauts during their 
interplanetary mission rely on radiation transport 
codes that accurately describe the interactions and 
propagation of these radiation fields throughout the 
bulk medium (ref. 1). HZETRN (ref. 2), which is 
a state-of-the-art galactic cosmic ray transport code 
developed at the Langley Research Center, contains 
detailed descriptions of physical processes such as 
energy loss, nuclear absorption, and fragmentation 
cross sections of projectile ions. However, the un- 
certainties of nuclear fragmentation cross sections 
existing in the code remain fairly large because of 
the lack of reliable experimental data (ref. 3). This 
condition is especially true for collisions of heavy 
ions on heavy ions for which the results of NUCFRG 
(ref. 4), based on a semiempirical abrasion- ablation 
model, are used as inputs to HZETRN. (Note that an 
improved version of a heavy ion fragmentation model 
will be available for use as inputs in the immediate 
future.) For collisions of heavy ions on protons, 
Rudstam’s semiempirical parameterization (ref. 5) 
has been used in the existing HZETRN. 

Rudstam and Metropolis et al. (refs. 5 and 6, 
respectively) were the pioneers who systematized 
high-energy cross section measurements into a use- 
ful analytical relationship describing nuclear reac- 
tions. These analytical relationships have been re- 


vised and improved by many new researchers as new 
cross-section data become available. The most com- 
prehensive set of semiempirical estimates of cross 
sections for a hydrogen target is perhaps due to 
Silberberg and Tsao (refs. 7 and 8), and to the most 
recent work of Webber, Kish, and Schrier (ref. 9). 
In the present study, the formulation of Silberberg 
and Tsao for a hydrogen target will replace that of 
Rudstam in HZETRN, and the effect of such a re- 
placement will be examined for estimates of galactic 
cosmic ray exposure using liquid hydrogen as a shield. 

Transport Methods 

As energetic ions traverse through bulk mat- 
ter, they lose energy through their interaction with 
atomic electrons along their trajectories. On occa- 
sion, they interact violently with nuclei of the matter 
and thus produce ion fragments or secondary nucle- 
ons moving in the forward direction and low-energy 
fragments of the struck target nucleus. The transport 
problem for the short-range target fragments can be 
easily solved in closed form in terms of collision den- 
sity (ref. 10) and treated separately. Hence, the 
projectile-ion fragment or secondary nucleon trans- 
port is the remaining problem of interest. In the 
previous work (ref. 10), the projectile-ion fragments 
and secondary nucleons were treated as if all went 
straight ahead (ref. 11). The straight-ahead approx- 
imation is found to be quite accurate for the nearly 
isotropic cosmic ray fluence (ref. 10). 


Because of the long range of the coulomb force and the large percentage of the material volume being 
occupied by electrons, the electron interactions can, to a good approximation, be treated as a continuous 
slowing-down process over any finite path length. Although the energy lost by an ion over some fixed path 
length fluctuates about a mean value, this fluctuation amounts to no more than a few percent (refs. 12 15) 
and is of no importance in the study of space radiation. In the following, a continuous slowing-down theory 
will be assumed, and the relevant quantity is the average loss per unit of path length. 

With the straight-ahead approximation and the target secondary fragments neglected (refs. 10 and 11), the 
transport equation may be written as 


~d__ 0_ 

dx Vj dE 


S(E) + a j(E) 


<t>j(x,E) = £ J™f jk (E,E') 4> k {x ,E') dE' 


( 1 ) 


where 4>j{x, E) is the flux of ions of type j with atomic mass Aj and charge Zj at x moving along the ar-axis at 
energy E (in units of MeV/amu), aj is the corresponding macroscopic nuclear absorption cross section, S(E) is 
the stopping power (continuous slowing-down approximation) of the protons, fjk{E, E r ) is a differential energy 
cross section for production of ion j in collision by ion k y and Vj is the range scaling parameter that is defined 
as 



(2) 


The solution to equation (1) is found to be subject to the boundary condition at x = 0 (that is, <5> 7 (0, E)), 
which is the incident beam spectrum. 

By transforming the heavy ion transport equation to an integral along the characteristic curve of that 
particular ion (ref. 10) and using the perturbation theory (ref. 16), the solution to equation (1) is given as a 
stepping procedure with step size h in the ^-direction (ref. 17). Thus, 

ipj(x + h, r) as exp[-<j(r, /i)] ipj(x, r + Vjh) 

+ Y1 f [ cx P[-<j ( r > z ) - (k( r \ h - z)\ 7 ) k {r + VjZ, r') 

^ JO Jr 

x ip k [x, r + u k {h - z)\ dr dz (3) 


where 

%!>j{x,Tj) = S{E ) 4>j(x, E) 

Tj is the residual range of ion j given by 



dE' 

sJW) 


and the exponential is the integrating factor with 


< j{r, t) = cr.j{r + Vjt') dt' 

Currently, we assume for Zj > 1 and k > j that 

Jjk( r > r ') — fjk( r ') *(r - r) 


( 4 ) 
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Using equation (4), equation (3) now becomes 

r}>j{x + h, r ) « cxp[-Cj(r, h)] ipj{x, r + Vjh) 

+ E / dz ex P[~Cj( r > z ) ~ C k( r '> h ~ z )] a .jk( r> ) 

k Jo 

xip k [x,r' + i>k{h- z)\ ( 5 ) 


with r' = r + ujz. Equation (5) is further approximated as 


ipj{x + h, r ) % exp[-Cj(r, h)] ipj{x, r + vjh) 


E j dz cxp[-C ; (r, 2 ) - Ct ( r > fl ~ z )\ a jk( r ) i’k{ x i r + + "feC* 1 - 2 )1 


k 

exp[— <Tj{r) h] tl>j(x,r + i/jh) 

f cxp[— <Tj(r) h] - cxp[-( 7 fc(r) /t] ) , 

°k(T)-Vj{r) J ’ 


+E^( r ) 

A: t 

+ 0[(^fc — ^j)^] 


r + ja/i) 


( 6 ) 


Equation (6) is the stepping formalism with energy-dependent cross sections for k > A He. The corresponding 
stepping formalism for nucleons has been discussed in detail in references 8 19. These stepping formalisms are 
then used to march the solution from the surface boundary to the desired shield thickness. 


Nuclear Fragmentation Data Base 

Even though the accuracy of the experimental 
data may improve for specific reactions, a reason- 
able means of representing data in computational 
procedures for cosmic ray transport calculation is 
still a challenge. For HZETRN and other radia- 
tion transport codes developed at Langley (refs. 18 
and 20), a point representation of the data is avoided 
because large multidimensional arrays will eventu- 
ally rival computer storage. Instead, various semi- 
empirical methods suitable for certain target or frag- 
ment groups arc put together in generating the data 
base. These scmiempirical methods, which are built 
on available experimental data and some theoretical 
base describing approximately the systematic varia- 
tion of reaction cross sections, offer the possibility 
of implementing any additional necessary correction 
factors or adjusting some existing parameters. 

When a nucleus is collided by high-energy nucle- 
ons, some individual nuclear constituents are ejected 
by direct knockout (ref. 21). The remaining nuclear 
structure is left in an excited state which seeks an 
equilibrium minimum-energy configuration through 
particle emission (ref. 5). This state is the basis of 
Rudstam’s formulation for the systematics of spal- 
lation products produced in such collisions. In his 
formalism, the distribution of resultant isotopes that 
arc related to the statistical nature of the evaporation 
process (ref. 22) is assumed to bo Gaussian centered 


at the nuclear stability line, and the total change 
in nuclear mass and the dependence on the incident 
projectile energy are treated empirically. Following 
Rudstam’s work, Silbcrberg and Tsao (refs. 7 and 8) 
later added more corrective factors to the formalism 
as more experimental data became available. In ad- 
dition, they added a scaling factor that relates the 
fragmentation cross section produced on a hydrogen 
target to that on the heavier target. 

Samples of fragmentation cross section aj k used 
as inputs to IIZETRN are shown as a function of 
projectile-ion charge Z k and fragment-ion charge Zj 
in figures 1 and 2 for a liquid hydrogen target and 
in figures 3 and 4 for a water target at various en- 
ergies of the projectiles. Figures 1 and 3 reflect the 
results produced by using Rudstam’s formalism for 
the fragments with Zj > 2, whereas figures 2 and 4 
reflect those by using the Silbcrbcrg-Tsao formal- 
ism. For nucleon and helium fragments, the cross 
sections were obtained from the results of Bcrtini 
et al. (refs. 23 and 24). (Note that the cross section 
shown in the figures has been reduced, for conve- 
nience, by a factor of 10 for nucleon fragments and a 
factor of 4 for helium fragments.) Because the scal- 
ing factor in the formalism of Silbcrberg and Tsao 
does not give adequate results for a target heavier 
than hydrogen (ref. 25), the scmiempirical abrasion- 
ablation model (ref. 4) is used in obtaining the cross 
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sections contributed by the oxygen ions in the water 
target. 

For all the cross sections shown in figures 1-4 
at three separate projectile energies (150, 600, and 
2400 MeV/amu), the peak of the fragments spectra 
(at constant Z *.) is generally higher as predicted 
by Rudstain than that predicted by Silberberg and 
Tsao. Conversely, the spread from the peak to 
the fragments with lower Z is lower for Rudstam, 
whereas the zig-zag shape over even and odd Z 
charges is more pronounced f«;r Silberberg and Tsao. 
At energies below 150 MeV/amu. '.ome Jesuits (not 
shown) given by Rudstam ’s formalism are erroneous; 
thus, for the transport calculations presented herein, 
the values of cross section at 150 MeV/amu have been 
extrapolated to below 150 McV for projectile ions 
with Z > 20, as has been done in reference 17. 

Results 

By using Rudstam’s formalism, several studies 
wore made in the past that involved galactic cosmic 
ray (GCR) transport calculations through a shield 
that contained hydrogen atoms, such as water (tis- 
sue equivalent) (refs. 26 and 27) or liquid hydrogen 
(ref. 28). In these calculations, the cross sections 
were assumed to be energy independent and set equal 
to the asymptotic values at high energy. Recently, 
improvements were made to IIZETRN by removing 
the energy- independent assumption; and the effects 
of such changes to the existing results were exam- 
ined in terms of LET spectra that showed a substan- 
tial enhancement of high- LET components (ref. 17). 
In the present study, Rudstam’s formalism was re- 
placed by that of Silberberg and Tsao, and the effect 
of that replacement was then examined by comparing 
the GCR exposure levels behind the liquid hydrogen 
shield at the solar minimum condition given by the 
CREME model of reference 29. 

The high-LET radiation components are usually 
degraded to lower LET components as a result of 
nuclear interactions between projectile and target 
nuclei, and such processes become more significant 
as the particles penetrate farther into the shield 
medium. This degradation is illustrated in figure 5 
in which the annual differential dose and dose equiv- 
alent arc plotted as a function of LET (or L) for 
2-, 5 and 15 g/cm 2 thicknesses of the shield. The 
spikes seen in the figures correspond to the zero slope 
of stopping power ( dS/dE = 0) at minimum ioniza- 
tion (ref. 30) of each ion, with protons starting at the 
lowest LET followed by increasing Z for the increas- 
ing L. Because the LET coordinate is plotted on a 
logarithm scale, the differential dose L dN/dL is con- 
verted to L 2 dN/dL (where N = f <p dE ) so that the 


area under the curve is linearly proportional to the 
total dose. Similarly, the differential dose equivalent 
is plotted as Ql? dN/dL , where the quality Q is a 
function of L. The magnification in dose equivalent 
at the high-LET region is a result of the ICRP 60 
quality factor (ref. 31). 

The differences between the spectra for the frag- 
mentation cross section by Rudstam’s formalism 
(fig. 1) and those by Silberberg and Tsao (fig. 2) 
as discussed earlier arc reflected in the calculated 
LET spectra for the transported flux through the 
liquid hydrogen shield. In figure 6(a), the ratio of 
the LET spectra for the differential flux calculated 
by using Silberberg and Tsao relative to that by us- 
ing Rudstam is displayed for the three thicknesses 
of the shield. The Silberberg- Tsao model produces 
more lower Z fragments and fewer higher Z frag- 
ments, and therefore more lower LET components 
and fewer higher LET components according to the 
spikes identified for each ion in figure 5. As a re- 
sult, the ratio (fig. 6(a)) is higher at LET values near 
the 100-MeV/cm region and lower at higher values of 
LET. These deviations from unity become more pro- 
nounced as the shield thickness increases. A similar 
comparison of integral flux LET spectra is shown in 
figure 6(b) in which the ratio at the lowest LET is al- 
most identical to unity; this indicates that the total 
flux does not change appreciably because of differ- 
ences in the fragmentation model. 

In the previous study (ref. 17), the removal of 
the energy-independent assumption w r as important 
in the risk assessment of GCR exposure. Using 
the liquid hydrogen shield as an example, the LET 
components above 100 keV//xm (or 1000 McV/cm), 
which could contribute some biological- risk orders 
of magnitude higher than the lower components, 
w r ere showm to increase by 40 percent at a 15-g/cm 2 
thickness because of the removal of the assumption. 
This conclusion was based on the use of Rudstam’s 
formalism. Similar results w^ere obtained in figure 7 
with the data base given by Silberberg and Tsao. 
Note in figure 5 that the region near 10 000 ke V / fim 
was not as critical because of the diminishing flux in 
this region. 

A comparison is made with the transport calcu- 
lations using the energy-dependent Silberberg- Tsao 
data base and the energy-independent data base from 
Rudstam to serve as a reference for the earlier GCR 
exposure studies made whth the old data base and 
the assumption of an energy-independent cross sec- 
tion which is taken to be the asymptotic value at 
2 GcV/amu. The ratio of the tw r o calculations in 
LET spectra, for both differential and integral flux 
shown in figures 8(a) and 8(b), respectively, displays 
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almost an oscillatory behavior over a wide range of 
LET. Thus, some of the differences are probably can- 
celled out. One arrives at the interesting observa- 
tion that an energy-independent calculation using 
the Rudstam cross sections yields results similar to 
those of an energy-dependent calculation using the 
Silberberg-Tsao values. 

Concluding Remarks 

Nuclear fragmentation cross sections of Silber- 
berg and Tsao that are more accurate for a hy- 
drogen target have been placed in the data base 
for calculations of galactic cosmic ray (GCR) trans- 
port. When compared with the old data base of 
Rudstam, the Silberberg-Tsao model produces fewer 
higher charge fragments but more lower charge frag- 
ments. Sample calculations of GCR transport with 
a liquid hydrogen shield reflect such differences of 
cross sections in that the transported flux based on 
the Silberberg-Tsao model contains fewer high-LET 
(linear energy transfer) components but more lower 
LET components. This disparity deepens as the 
shield thickness increases. When the Silberberg-Tsao 
cross sections are assumed to be energy indepen- 
dent in the data base, the comparative results in- 
dicate an underestimation of the LET components 
above 100 keV//im by approximately 40 percent for a 
15-g/cm 2 thickness of liquid hydrogen caused by the 
assumption. Moreover, the present study found that 
an energy-independent calculation would be best ac- 
complished by using Rudstam’s cross sections as done 
in the past. 

NASA Langley Research Center 
Hampton, VA 23681-0001 
June 8, 1993 
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(c) Projectile energy = 2400 MeV/amu. 

Figure 1. Semiempirical, differential fragmentation cross section for liquid hydrogen target as function of 
projectile- ion charge and fragment-ion charge Zj according to Rudstam’s formalism for fragments heavier 

than helium. Cross sections shown for nucleon and helium fragments have been reduced by factors of 10 
and 4, respectively. (1 cm 2 /g converts to 1.67 barns.) 
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Gjb cm2 /g 




Oj k , cm2/g 


(c) Projectile energy = 2400 MeV/amu. 

Figure 2. Semiempirical, differential fragmentation cross section for liquid hydrogen target as function of 
projectile-ion charge and fragment-ion charge according to Silberberg-Tsao formalism for fragments heavier 
than helium. Cross sections shown for nucleon and helium fragments have been reduced by factors of 10 
and 4, respectively. (1 cm 2 /g converts to 1.67 barns.) 
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(a) Projectile energy = 150 MeV/amu. 


Cjb cm2 /g 



(b) Projectile energy = 600 MeV/amu. 



(c) Projectile energy = 2400 MeV/amu. 

Figure 3 Semiempirical, differential fragmentation cross section for water target as function of projectile- 
ion charge and fragment-ion charge according to Rudstam’s formalism for fragments heavier than helium. 
Cross sections shown for nucleon and helium targets have been reduced by factors of 10 and 4, respectively. 
(1 cin 2 /g converts to 29,9 barns.) 
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ffjb cm2 /g 


(a) Projectile energy = 150 MeV/amu. 


°jb cm 2 /g 


(b) Projectile energy = 600 McV/amu. 




Figure 4. Scmiempirical, differential fragmentation cross section for water target as function of projectile- 
ion charge and fragment-ion charge according to Silberberg-Tsao formalism for fragments heavier than 
helium. Cross sections shown for nucleon and helium fragments have been reduced by factors of 10 and 4, 
respectively. (1 cm 2 /g converts to 29.9 barns.) 
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L, MeV/cm 


(a) 2-g/cm 2 thickness. 



(b) 5-g/cm 2 thickness. 

Figure 5. Differentia] LET spectra for annual dose equivalent behind various thicknesses of liquid hydrogen 
shield exposed to GCR at solar minimum. 
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(c) 15-g/cm 2 thickness. 
Figure 5. Concluded. 




L, MeV/cm 

(a) Differential LET spectra. 



L, MeV/cm 

(b) Integral LET spectra. 

Figure 6. Comparison of LET spectra for transmitted GCR flux calculated by using Silberbcrg-Tsao energy- 
dependent cross sections relative to the flux by Rudstam through several thicknesses of liquid hydrogen 
shield at solar minimum. 
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(a) Differential LET spectra. 



(b) Integral LET spectra. 

Figure 7. Comparison of LET spectra for transmitted GCR flux calculated by using Silberberg-Tsao energy- 
dependent cross sections relative to the flux by energy-independent cross sections through several thicknesses 
of liquid hydrogen shield at solar minimum. 
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(a) Differential LET spectra. 



uJ 

105 


Figure 8. Comparison of LET spectra for transmitted GCR flux calculated by using Silberberg-Tsao energy- 
dependent cross sections relative to the flux by Rudstam energy-independent cross sections through several 
thicknesses of liquid hydrogen shield at solar minimum. 
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